SSHPU 

f 


NASA Contractor Report 185197 


Viscous Three-Dimensional 
Analyses for Nozzles for 
Hypersonic Propulsion 


G.J. Harloff, D.R. Reddy, 
and H.T. Lai 

Sverdrup Technology, Inc. 

NASA Lewis Research Center Group 
Cleveland, Ohio 


January 1990 


Prepared for 
Lewis Research Center 

Under Contracts NAS3-24105 and NAS3-25266 




National Aeronautics and 


Space Administration 


(NASA— CR— 135 197) VISCOUS 
ANALYSES FOR NOZZLES FOR 
PROPULSION Final Report 
Technology) 25 p 


THREE-OIMFNS IONAL 

hypersonic 

( Sverdrup 

CSCL 21E 


G3/07 


N90-17636 


Unci as 
0264834 


E-5267 


VISCOUS THREE-DIMENSIONAL ANALYSES FOR NOZZLES FOR HYPERSONIC PROPULSION 


G.J. Harloff, D.R. Reddy, and H.T. Lai 
Sverdrup Technology Inc. 

NASA Lewis Research Center Group 
Cleveland, Ohio 44135 


SUMMARY 

A Navier-Stokes computer code has been validated using a number of two- and 
three-dimensional configurations for both laminar and turbulent flows. The validation 
data covers a range of freestream Mach numbers from 3 to 14, includes wall pressures, 
velocity profiles, and skin friction. Nozzle flowfields computed for a generic scram- 
jet nozzle from Mach 3 to 20, wall pressures, wall skin friction values, heat transfer 
values, and overall performance are presented. In addition, three-dimensional solu- 
tions obtained for two asymmetric, single expansion ramp nozzles at a pressure ratio 
of 10 consist of the internal expansion region in the converging/diverging sections 
and the external supersonic exhaust in a quiescent ambient environment. The fundamen- 
tal characteristics that have been captured successfully include expansion fans; Mach 
wave reflections; mixing layers; and nonsvmmetrical , multiple inviscid cell, super- 
sonic exhausts. Comparison with experimental data for wall pressure distributions at 
the center planes shows good agreement. 

NOMENCLATURE 

a sound speed 

Cf skin friction coefficient 

Ch heat transfer coefficient 

Cp specific heat 

C v velocity coefficient 

M Mach number 

NPR Nozzle Pressure Ratio 

n nozzle efficiency 

Pr Prandlt number 

P w local wall static pressure, Pw/P(o) 

Re Reynolds number 

T temperature 

X streamwise distance 

Y ratio of specific heats • 

Subscripts : 
c combustor exit 

e nozzle exit 

1 laminar 

n unit normal 

o freestream condition 

s static condition 

t total condition 

t turbulent 

w wall 

vx viscosity 

p density 



INTRODUCTION 


Accurate evaluation of nozzle performance is essential to hypersonic propulsion 
because that performance is highly sensitive to net thrust. Highly integrated 
fuselage/nozzle configurations may experience complex interaction of shocks, turbulent 
mixing, differing levels of under/over expansion, and possible boundary layer separa- 
tion. To simulate these complex interactions, the full Navier-Stokes equations need 
to be solved and the appropriate computer code carefully calibrated before it can be 
applied to the hypersonic nozzle flow field. The goal of the current study is to vali- 
date the PARC3D code over a range of flowfield conditions and to examine the flow phys- 
ics and overall nozzle performance over a wide range of flight conditions. The PARC3D 
code solves the full, three-dimensional Reynolds averaged Navier-Stokes equations in 
strong conservation form with the Beam and Warming approximate factorization algorithm; 
the implicit scheme uses central differencing for a curvilinear set of coordinates. 

The code was originally developed as AIR3D by Pulliam and Steger; 1 Pulliam later added 
Jamison's artificial dissipation and called it ARC3D . 2 Cooper 2 then adapted the ARC3D 
code for internal propulsion and renamed the code PARC3D . 


Code Validation 

To gain confidence with the PARC 3D code simulating the interaction phenomena 
characteristic of hypersonic nozzle flow, experimental data were selected for code val- 
idation as follows: a two-dimensional shock boundary layer interaction at Mach 14 
(Holden 4 ), and a Haack-Adams body (Allen 5 ) with a favorable pressure gradient at 
Mach 3. Previous configurations examined at NASA Lewis by Reddy and Harloff 6 include: 
a two- and a three-dimensional hypersonic nozzle at Mq = 5, for which comparisons were 
made with computations of Spradley 7 and others; and a three-dimensional corner flow 
configuration at Mq = 2.94 of Oskam et al. 8 The Reddy and Harloff 6 study provided 
confidence in the PARC codes for a variety of complex shapes and high speeds. 


Two-Dimensional Shock-Boundary Layer Interaction at Mach 14 

The two-dimensional experimental data 4 at Mach 14 are used as the first test 
case. The flat plate length was 1.44 ft, the ramp angle was 18°, and the ramp length 
was 1.14 ft. The Reynolds number per foot, based on total conditions, was 32xl0 6 and 
the computational grid 399 x 99. The predicted surface pressure coefficient and skin 
friction are compared to the test data (Figs. 1 and 2); good agreement for both is ev- 
ident. The predicted negative skin friction for the separated zone is also in good 
agreement with experiment. As the boundary layer in the experiment was laminar, this 
study validates the code at hypersonic conditions for laminar flow. 


Haack-Adams Body at Mach 3, Turbulent 

The axisymmetric Haack-Adams body, 5 was 3 ft long, had a length to diameter ratio 
of 10, and a base to maximum diameter ratio of 0.532. The freestream Mach number was 
2.96 and the Reynolds number based on static conditions O.lxlO 6 cm. Detailed velocity 
measurements were made at 7 axial stations. Comparison of measured and computed veloc- 
ity profiles at stations 2 and 7 (Figs. 3 and 4) are in excellent agreement with each 
other. Axial stations shown in Figs. 3 and 4 were at 0.28 and 0.97 of the body length, 
respectively. The streamwise pressure gradient is favorable for this configuration, 
which is similar to the nozzle flow environment. The computed Cf (Fig. 5) shows ex- 
cellent agreement with experimental data. This study validates the 3aldwin-Lomax tur- 
bulence model used in the PARC3D code for attached axisymmetric flow at a single Mach 
number . 


NOZZLE FLOWFIELD CALCULATIONS 


Two-Dimensional Flowfield Calculations, Mach 5, Laminar 

A generic three-dimensional nozzle 7 is shown schematically in Fig. 6. Its length 
was six entrance nozzle heights, the upper wall has a 20° slope, the lower splitter 
length was three. The lower wall was horizontal up to one, where the wall expands at 
a 6° angle. The flow was assumed to be laminar. The grid used was 199 x 99. 


The nozzle entrance flow was assumed to be uniform, and the ratio of specific 
heats 1.27 for both streams, the nozzle entrance velocity 1610 m/sec., the static pres- 
sure 9206 N/m 2 , the static temperature 2311 °K, and the Mach number 1.657. The corre- 
sponding freestream values were: 1765 m/s ec , 506 N/m 2 , 261 °K, and 5.0, respectively. 
Upper wall static pressures and values from the GI.M code, the Seagull code, 2 and the 
Method of Characteristics (MOC) 9 are shown in Fig. 7. Pressures near the nozzle en- 
trance are in agreement with the GIM results, but are higher in the aft region. The 
inviscid results of Seagull and MOC lie between the current viscous solution and the 
GIM solution in this region. The predicted flowfield is presented in terms of Mach 
number in Fig. 8. The contact surface between the nozzle flow and the freestream de- 
flects downward, as is evident in the total pressure contour plot and the velocity vec- 
tors. 9 Three-dimensional calculations were also performed for this nozzle. 6 
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Two-Dimensional Flowf ield , Mach 3 to 20 , Turbulent 

Figure 9 shows the nozzle geometry used in this study. The initial divergence 
angle is about 37°. 


Freestream and Combustor Exit Properties 

The ratio of specific heats, y, and the gas constant of the combustor exit are 
considered constant throughout the entire flowf ield (Table I). Because the static noz- 
zle entrance Reynolds numbers, based on nozzle entrance height, are such that turbu- 
lent flow is expected for the entire range of Mach numbers studied, turbulent flow is 
assumed in the whole flowfield. 


GRID 


The grid used was 199 x 129, uniform in the x-direction (Fig. 10). Grids were 
clustered by being stretched by hyperbolic functions near the upper and lower walls to 
resolve the boundary layers, as well as downstream of the cowl to resolve the shear 
layer. The grid was not adapted to the shear layer location, although this has been 
accomplished by the authors in previous unpublished studies, which indicated that the 
present grid was sufficient to capture the flowfield, 

Boundary conditions are uniform inflow and extrapolated outflow properties for 
supersonic flow. The supersonic/hypersonic inflow boundary above the exit was assumed 
to be uniform. For subsonic outflow, the pressure is imposed by the sublayer approxi- 
mation. These boundary conditions are ideal and more realistic inflow profiles can be 
accounted for. The wall temperature is assumed to be constant at 1750 °R. 


Mach 3 Flowfield 

Velocity vectors for Mq = 3 (Fig. 11) show that the flow is overexpanded and a 
shock propagates from the cowl lip to the upper nozzle wall. The exhaust flows along 
the upper nozzle wall similar to a wall jet. A slight separation occurs on the upper 
wall ahead of the shock (discussed below in terms of Cp) . The Mach number contours 
are shown in Fig. 12. Due to the overexpanded flow at the find of the lower cowl, the 
shear layer turns upward, which causes a shock to form above it. 

Total pressure contours shown in Fig, 13 illustrate the gradient in nozzle and 
freestream total pressure across the shear layer. 

The divergence of velocity provides a scalar measure of density gradients such 
that compressive and expansive regions of the flow can be visualized as dark and light 
regions respectively (Fig. 14). The deflection angle of the shock from the cowl end 
to the upper nozzle wall is about 39° and the Mach number in front of the oblique 
shock is about 3.5. After the shock reflects off the boundary layer, it propagates in 
the downstream direction. The boundary layer separates on the upper nozzle wall ahead 
of the incident shock (further illustrated with the skin friction prediction below) . 
Details of this region show a small region of reverse flow. As the shock interacts 
with the wall boundary layer, it initially reflects at a sharp angle. As the flow 
moves over the surface of the separation bubble, it is turned less, as shown by the de- 
crease in magnitude of the compression. At the aft end of the separation bubble, the 
flow is again turned at a sharp angle and coalesces into the reflected shock down- 
stream of the separation. Below the cowl, a shock emanating from the leading edge of 
the cowl is observed, and is caused by flow deflection from the boundary layer along 
the lower cowl wall. The expanding region of the flow downstream of the cowl is also 
observed . 

Wall pressures as a function of X for the three nozzle wall surfaces 
(Fig. 15(a)) illustrate that the wall pressure initially decreases with increasing X 
for both upper and lower nozzle walls. The interior nozzle shock intersects the upper 
nozzle wall at about 2.8 ft, as evidenced by the marked pressure rise there. Because 
the pressure beneath the cowl is greater than the pressure in the nozzle, the shock 
forms in the nozzle. 

Wall skin friction values are shown in Fig. 15(b) for the three surfaces. The 
skin friction on the upper nozzle wall becomes negative at about 2,8 ft, indicating a 
small region of reverse flow. The skin friction of the upper cowl wall is slightly 
higher than that of both the upper and lower walls. The heat transfer for the upper 
and lower nozzle walls is similar over the cowl's length (Fig. 15(c)) while the rise 
in the heat transfer coefficient is observed at X = 2.8 ft, where the boundary layer 
flow separates. In the nozzle the heat transfer is positive, indicating hot gas and a 
cold wall, but below the cowl the heat transfer coefficient is negative, indicating 
the reverse. 

A large scale separation on the upper nozzle surface is predicted for laminar 
flow (Fig. 16), whereas the turbulent computations show a small separation. 
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Mach 6 Flowfield 


Figure 17 shows the velocity vectors for Mq = 6 . The flow, slightly overex- 
panded at the cowl exit, expands behind the cowl by means of an expansion fan. The 
wake is deflected upward, as indicated by the velocity vector angles and the position 
of the velocity shear layer. No reverse flow at the upper boundary is indicated. 

Mach number contours for the flowfield (Fig. 18) also indicate the location of the 
shear layer while Fig. 19 shows its upward deflection as well as total pressure con- 
tours . 

Wall pressures (Fig. 20(a)) show that the interior wall pressure is overexpanded. 
The wall skin friction and heat transfer values are shown in Figs. 20(b) and (c) 
respectively. 


Mach 10 Flowfield 

The exhaust flow is underexpanded at Mach 10 as indicated by the downward deflec- 
tion of the nozzle flow (Fig. 21). There are no indications of shocks inside the noz- 
zle, but a shock propagates from the lower cowl lip downward due to the underexpansion 
of the flow (Fig. 22). The total pressure contours at Mq = 10 (Fig. 23) show that 
the flow is underexpanded and that the shear layer is deflected downward. The shock 
deflection angle is about 18°. Pressure is constant on the lower wall (Fig. 24(a)), 

Cf and Cft wall values are also shown in Fig. 24. 


Mach 20 Flowfield 

The velocity vectors are deflected downward aft of the cowl due to under expans ion 
(Fig. 25) which causes a shock to propagate from the lower cowl downward. The lower 
cowl leading edge shock merges with the trailing edge shock about 26 percent of a cowl 
length downstream of the cowl. 

The flow along the curved upper nozzle surface turns and compresses, like flow 
along a curved ramp, and a curved shock forms inside the nozzle. An upwardly deflected 
shock behind the end of the upper nozzle wall is evident (Figs. 25 and 26). The Mach 
number contours are shown in Fig. 26. Total pressure contours (Fig. 27) show the down- 
ward deflection of the shear layer. The shock is detached from the lower side of the 

cowl due to the boundary layer there. The shock from the leading edge of the cowl 
merges with the shock from the aft end about 18 percent of the length downstream of 
the cowl. An interior shock along the upper wall turns the flow towards the horizon- 
tal axis. This interior shock is typical at overspeed conditions. The nozzle flow is 
underexpanded (Fig. 28(a)); the upper wall pressure does not follow the lower nozzle 
wall pressure but instead increases the full length of the nozzle, even inside the 
cowl/ The Cf and C p values are shown in Fig. 28(b) and (c) respectively. 

A perfect gas model has been assumed in these calculations, and the constant ratio 

of specific heats (y) has been varied with flight Mach number. This assumption may 
not be as valid at high Mach numbers as at low Mach numbers and should be further 
investigated . 


Nozzle Performance 

The adiabatic nozzle efficiency, N, is the ratio of the nonisentropic to isentro- 
pic expansion enthalpies between the combustor exit, c, and the nozzle exit, e. The 
velocity coefficient, C v , is the ratio of the actual to the isentropic velocity and is 
the square root of the efficiency. The velocity coefficient increases with Mach 
number up to Mq of 10 but decreases at Mq of 20. For example, C v is 0.93 at Mq 
of 3, 0.98 at M 0 of 6, 1.0 at M 0 of 10 and 0.94 at M 0 of 20. The reported C v 
values peak at Mq of 10 because the nozzle design Mach number is close to 10 and 
the Mq 20 is an overspeed condition for this nozzle. The frozen gas assumption has 
been successfully used for ramjet, rocket, and hypersonic wind tunnel nozzle flows. 
Because these values might be different with chemistry accounted for, the above C v 
values should be considered as relative. For example. Ref. 11 has reported that 
finite-rate chemical recombination of dissociated species affects scramjet nozzle 
force coefficients, at Mach 20, and that the nozzle flow analyzed was in vibrational 
equilibrium . 


Three-Dimensional, Single Expansion Ramp Nozzle, NPR = 10, Laminar 

Three-dimensional solutions of a single expansion ramp nozzle are performed to 
simulate the nonaxisymmetr ic nozzle flowfield in both internal/external expansion re- 
gions and the exhaust plume in a quiescent ambient environment. Two different configu- 
rations of the nozzle at NPR - 10 are examined. They were selected from the single 
expansion ramp nozzles experimentally investigated by Re and Leavitt 11 to analyze the 
effects of various geometrical parameters and pressure ratios on the nozzles 1 
static performance. 

The present converging/diverging nozzle has a long upper surface functioning as an 
external expansion ramp, with a rectangular cross-section in the internal nozzle. The 
edge of the flat sidewall is highly skewed at the nozzle exit to connect the upper and 
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lower: surface lips. The flow generated by large streamwise pressure differences ex- 
pands from near reservoir conditions through the nozzle and exhausts supersonically 
into the quiescent air. The resulting flowfield is characterized by flow expansion 
inside the nozzle and by exhaust plume interaction with the ambient environment. In 
particular, a strongly interactive flow structure containing a normal shock, expansion/ 
compression wave reflection, and separation can occur along the external expansion 
ramp of the upper surface. Around this section, the free shear layer emerging from 
the" lip of the shorter lower surface acts as an artificial nozzle wall, and its trajec- 
tory shape can alter the nature of the flow, depending on the ambient conditions below 
it. The wave structure initiated by the shear layer may include an oblique shock or 
an expansion fan emanating from the lower nozzle lip. This shear layer can also de- 
flect. the flow upwards at a high angle to interact with the upper boundary layer. 

A similar pattern exists in the spanwise direction, where the sidewall creates a 
vertical free shear surface starting at the skewed nozzle exit. At high Reynolds 
number, the shear layer behind a sharp or blunt trailing edge may not be stable; how- 
ever high Reynolds number interaction of the free shear layer is even more unstable in 
a quiescent external environment than in a subsonic or supersonic external stream. 
Consequently, a steady state solution may not be possible for certain flows with a 
quiescent external stream. The three-dimensionality of the flowfield, its unsteadi- 
ness, and stability considerations are some of the complications of the present prob- 
lem, although only steady-state solutions are computed. 

Converged solutions are compared with the experimental measurements for the 
pressure distributions on the upper and lower nozzle surfaces. The configuration of 
the nozzle (Fig. 29) has an upper two-dimensional flap and a shorter, lower two- 
dimensional flap, both flat in the spanwise direction. Part of the upper flap extend- 
ing from the throat area functions as an external expansion ramp. The edge of the ver- 
tical sidewall is highly inclined. This arrangement gives a rectangular cross section 
parallel to the exit plane aligned with the tilted edge. Intersections of the tilted 
edge and the upper and lower nozzle walls normally occur upstream of the ends of the 
wall surfaces . 

Two nozzle configurations are examined (Fig. 30) with identical lower flaps and 
sidewalls and a nozzle aspect ratio of width to height of 4, with the height measured 
at the throat. Note that the flow expands to a lower pressure in case 2 because the 
external expansion ramp is longer. The flat sidewall is assumed to have a very thin 
thickness (0.0058 cm). 

All boundary conditions are treated explicitly. Since a plane of symmetry exists 
in the spanwise z-direction, only half of the nozzle flowfield is computed. Symmetry 
conditions are specified on this centerplane. No-slip velocity and an adiabatic wall 
temperature are imposed on the nozzle surfaces. In the farfield, still air conditions 
at 101.3 kPa and 300 °K, are fixed at the upper and lower boundaries in the vertical 
y-direction sufficiently far away from the nozzle; similar boundary conditions are 
imposed at the farfield in the spanwise direction. In the streamwise x-direction, a 
stagnation pressure of 1013 kPa and a stagnation temperature of 300 °K are specified 
at the nozzle entrance. Other variables at this location are computed using the isen- 
tropic relationships and a characteristic-like condition extrapolated from the inte- 
rior, At the outflow boundary, streamwise flux gradients are assumed to be negligible. 

For normalization, stagnation quantities at the nozzle entrance are taken as the 
reference variables. The nozzle throat height and the speed of sound are the refer- 
ence length and velocity, respectively. The Reynolds number is computed based on 
these values. The reference temperature, pressure and length are taken to be 300 °K, 
1013 kPa, and 2.54 cm, respectively. The corresponding Reynolds number is 5.618 mil- 
lion., 


Grid 


The simulated flow includes the nozzle interior and exterior, which contains the 
downstream exhaust. Since the nozzle has a symmetrical plane at the middle in the 
spanwise direction, the spanwise regions include the interior half of the nozzle and 
the flow exterior to the sidewall. The complete domain' is then divided correspondingly 
into subsections to simplify grid generation using an algebraic technique. To resolve 
the viscous layers, hyperbolic tangent functions are used to cluster grid points in 
regions near the walls. Patching these separate segments produces a final sheared grid 
in Cartesian coordinates. 

Figure 31(a) illustrates a view of the nonuniform grid in the x-, y-plane along 
the streamwise direction. The grid clusters near the wall surfaces to the exhaust 
region behind the inclined nozzle exit. Figure 31(b) presents a detailed view of the 
same grid showing a high degree of grid nonorthogonality, particularly next to the 
exit. The skewed grid lines result since one of the vertical curvilinear coordinates 
aligns with the sidewall tilted edge to simplify the application of the boundary condi- 
tions. Spanwise views of the grid distribution in the y- z plane are shown at two dif- 
ferent axial locations inside the nozzle (Fig. 32(a)) and downstream in the exhaust 
plume, (Fig. 32(b)). The vertical coordinates cluster about the interior and exterior 
surfaces of the sidewall, and are not on a plane with constant x-coordinate . 

The flow starts from the nozzle entrance, X = 0, where the total pressure is pre- 
scribed, expands through the converging/diverging and external ramp sections, and 
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exhausts into the downstream low/pressure region, which is bounded by quiescent bound- 
aries far from the nozzle. The farfield computational boundaries are approximated at 
distances of 25, 10, and 8 characteristic lengths in the x-, y- and z-direct ions , re- 
spectively. The solutions discussed below have been computed on a 135 by 120 by 50 
grid for both cases. Of these, 60 by 40 by 30 grid points are employed in the nozzle 
interior. The smallest grid is on the order of 0.001 near the wall regions for the v- 
and z-direct ions . This grid gives typically two to four subsonic points in the domi- 
nant supersonic viscous layers. For the x-direction, the grid is slightly refined 
near the nozzle exit and is on the order of 0.01. In the figures discussed below, i 
and 1 indicate the streamwise and spanwise grid indices respectively, while x and 
z represent the corresponding nondimensional coordinates. 


Case 1 : 

The converged result of the first case is presented with Mach number, density, 
and pressure contours at several streamwise and spanwise cross sections; contours are 
in equally spaced increments over the entire range of values unless otherwise noted. 
Mach number contours on the center plane (Fig. 33(a)) show a typical exhaust plume con- 
taining a steady-state, multiple inviscid cell pattern in the supersonic exhaust 
core. The flow is sonic at the throat and highly supersonic at the exit. The first 
inviscid cell appears very irregular and is bounded by an expansion/compression wave 
system. Because of the inclined trajectory of the exhaust plume, the Mach waves ema- 
nate from the ends of the flap surfaces, where the upper and lower free shear layers, 
which depart from the nozzle walls at different angles, become the expansion and com- 
pression ramps. The lower shear layer is very diffusive and largely due to lack of 
grid resolution along this shear layer. The first cell has an average Mach number of 
2.7, while the Mach numbers near the cell centers of the next two cells are about 2.3 
and 2.0, respectively, showing a gradual decay of the core center velocities. The 
core cells reduce in size and numbers towards the sidewall, (Figs. 33(b) to (d)) whose 
effects become pronounced (Fig. 33(d)) with the appearance of the initial vertical 
shear layer. The flow inside the nozzle has very thin viscous layers next to the 
walls, particularly on the lower surface, and thus appears predominantly inviscid. 

The spanwise contour views are shown in Figs. 34(a) to (d) . Figure 34(a) 
presents Mach contours at an axial station in the nozzle area, where the clustered 
regions indicate the upper, lower, and sidewall surfaces. The contours in Fig. 34(b) 
are located on the external expansion ramp. The structure includes the lower shear, 
thin vertical shear, and upper boundary layers. In the exhaust region, these free 
shear layers become thicker with downstream distance (Figs. 34(c) to (d)) in which the 
plume boundaries are evident. The density contours at the center plane and on the 
plane next to the sidewall (Figs. 35(a) to (b)) show a pattern very much like Mach con- 
tours but out of phase. Density variation is also small across both shear layers.- 

This same variation occurs for the pressure (Fig. 36) which decays quickly to the 
ambient value. The contours (Figs. 35(a) and 36) show a weak wave system in the 
exhaust core, particularly after the first inviscid cell. Figures 37(a) and (b) , and 
38 depict typical density and pressure contours in the spanwise direction, plotted with 
smaller increments for clarity. Note that large variations in density and pressure 
occur inside the nozzle, and the exhaust flow does not deviate much from the perfectly 
expanded condition. Finally, comparison made with the measured data (Figs. 39(a) to 
(b)) for the upper and lower surfaces in the center plane show good agreement for wall 
pressure distributions. 


Case 2: 

Solution of the second case is also presented in terms of the contours of the flow 
variables. Figures 40(a) to (d) depict the basic flow pattern of the supersonic 
exhaust into quiescent air. For this geometry, the upper plume boundary is well-defined 
by the very thin free shear layer which follows the upper nozzle contour, and remains 
nearly at this height downstream. The other shear layer is also deflected downwards 
and is diffusive. Although the flow expands to about the same maximum Mach number in- 
dicated in case 1, the supersonic inviscid core (Fig. 41(a)) is more extensive. As a 
result, the fluctuating decrease of the core center velocities becomes smaller. More- 
over, the flow is highly asymmetric along the exhaust center line. The sidewall 
(Fig. 41(d)) causes the initial development of the vertical shear layer, which curves 
irregularly towards the center plane as it progresses downstream since the spanwise 
velocity component is large and positive in the initial region of the vertical shear 
layer, but becomes small and negative downstream. The flow behind the inclined exit 
consists of two small, separate supersonic regions next to the exit plane. 

Spanwise views of Mach contours at various axial locations are given in 
Figs. 41(a) to (d) . The interior cross-sectional view (Fig. 41(a)) shows a thicker ex- 
pansion layer in the lower half of the nozzle. Figures (41(b) to (d)) depict the con- 
tours at the axial plane on the external ramp region and the planes in the exhaust 
region. Density and pressure contours in the streamwise planes are shown in 
Figs. 42(a) to (b), and 43. A clear indication of the compression/expansion wave re- 
flection can be observed at the first cell boundary Figs. 42(a) and 43 after which the 
wave structure then becomes very weak with downstream distance. Spanwise variations 
of density and pressure are illustrated in Figs. 44(a) to (b) and 45. Comparison 
between the numerical prediction and experiment (Figs. 46(a) to (b)) show very good 
agreement for the wall pressure distributions in the center plane. 
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SUMMARY OF RESULTS 


The PARC3D computer code has been validated using a number of two- and three- 
dimensional configurations with both laminar and turbulent flows. The validation data 
used included wall pressures, velocity profiles, and skin friction. The validation 
data cover Mach 3 to 14 with shock-boundary layer interaction and/or boundary layer 
separation . 

Two-dimensional nozzle flowfields were computed and overall performance for a re- 
presentative scramjet nozzle from Mach 3 to 20. Matching the external pressure at the 
cowl exit is necessary to avoid waves in the nozzle. At Mq of 3 and 6 the nozzle 
flow is overexpanded and shocks are inside the nozzle, while the shear layer is de- 
flected upward. The laminar Mg = 3 nozzle flowfield showed a massive separation on 
the upper ramp surface, whereas the turbulent solutions showed only a small separa- 
tion. At Mq of 10 and 20 the nozzle flow is underexpanded and shocks are outside the 
nozzle, while the shear layer is deflected downward. At Mq = 20, a shock forms inside 
the nozzle along the upper wall to deflect the flow from the radial to the horizontal 
direction . 

Three-dimensional solutions were obtained for two asymmetric, single expansion 
ramp nozzles at a pressure ratio of 10. The computed flow consists of the internal ex- 
pansion region in the converging/diverging sections and the external supersonic exhaust 
in a quiescent ambient environment. The fundamental characteristics of the flowfields 
for the prescribed flow conditions were predicted successfully and include expansion 
fans; Mach wave reflections; mixing layers; and nonsymmetr ical , multiple inviscid cell, 
supersonic exhausts. Comparison with experimental data for wall pressure distributions 
at the center planes shows good agreement. 
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TABLE I . - FREESTREAM AND COMBUSTOR EXIT PROPERTIES 


M 0 

T 0' 

Pq. 

Me 

T C ' 

Pc- 

Y 

R, 


R 


R 



ft-ibf/; 



psia 



psia 


Ibm °P 

3 

390 

1.7 

1.0 

4075 

11.9 

1.25 

62.7 

6 

393 

. 4 

1.5 

4220 

10.4 

1. 19 

63.4 

10 

409 

. 15 

1.9 

5403 

75 . 0 

1 . 16 

66 . 0 

20 

468 

. 03 

4 . 6 

5119 

25.9 

1 . 18 

88 . 0 


HOLDEN'S DATA 

COMPUTATION 

RANKINE-HU60NI0T 



HOLDEN'S DATA 
COMPUTATION 



FIGURE 1. - Cp ON ?8 DEGREE RAMP AT M - 1R.1. 


FIGURE 2. - Cf ON 18 DEGREE RAM 3 AT M 


EXPERIMENTAL 

COMPUTATION 



EXPERIMENTAL 

COMPUTATION 



500 1000 1500 2000 

VELOCITY/ Ff/SEC 

FIGURE 3. - HAACK - ADAMS BODY AT MACH 2.96. 
(STATION 2. X/L = 0.28). 


500 1000 1500 201 

VELOCITY, I T/SEC 

FIGURE <l. - HAACK - ADAMS BODY A1 MACH 
(STATION 7. X/L = 0.97>. 







FIGURE 9. - NOZZLE DEFINITION. FIGURE 10. - TYPICAL GRID FOR COMPUTATION 


UZERO 



FIGURE 11. - MACH 3, VELOCITY VECTORS. 


FIGURE 12. - MACH 3, MACH NUMBER CONTOURS 




FIGURE 13. - MACH 3, TOTAL PRESSURE CONTOURS. 



FIGURE 14. - MACH 3, DIVERGENCE OF VELOCITY. 
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FIGURE 15. - MACH 3, P w . C«, C h VERSUS X. 


„*vkkk* 

,kkk V V k k W+iyjWW 


UZERO 



FIGURE 17. - MACH 6, VELOCITY VECTORS. 
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B 0.71 
C 1.42 
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FIGURE 19. - MACH 6, TOTAL PRESSURE CONTOURS 






FIGURE 23. - MACH 10, TOTAL PRESSURE CONTOURS. 






FIGURE 25. - MACH 20 . VELOCITY VECTORS. 


FIGURE 28. - MACH 20. MACH NUMBER CONTOUR 


A 0.00 
B 1783.24 
C 3566.49 
D 5349.74 
E 7132.99 
F 8916.25 
G 10699.49 



FIGURE 27. - MACH 20. TOTAL PRESSURE CONTOURS. 
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— UPPE3 COWL WALL 
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AXIAL DISTANCE 

FIGURE 28. - MACH 20. P . C f , C h VERSUS X. 
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FIGURE 30. - COMPARISON OF T-JO NOZZLE CONFIGURATIONS 








(d) ! = 25 . z = 1.982. 

FIGURE 53. - STREAMWISE MACH NUMBER CONTOURS, CASE 1. 



(b) 1 = 29, z = 1.998. 

FIGURE 35. - STREAMWISE DENSITY CONTOURS, CASE 1. 


(a)i = 50, x = 4.485. (I>) I* 65. X s 5.372. 



(C) i - 75, X =• 6.221. (( )i = 85, X - 7.513. 

FIGURE 34. - SPANWISE MACH NUMBER CONTOURS, CASE 1 



FIGURE 36. - STREAMWISE PRESSURE CONTOURS, CASE 1,1-1 



FIGURE 37. - SPANWISE DENSITY CONTOURS, FIGURE 38. - SPANW1SE PRESSURE 



0 2 4 6 3 1C 

(b) LOWER WALL. 

FIGURE 39. - WALL PRESSURE DISTRIBUTIONS, CASE 1. 




FIGURE 43. - STREAMWISE PRESSURE CONTOURS, CASE 2, 1 
Z = O.O. 


(b) 1 = 29, Z = 1.998. 

FIGURE 42. - STREAMWISE DENSITY CONTOURS, CASE 2. 



(b) i = 90, x = 8.377. 
FIGURE 44. - SPANWISE 
DENSITY CONTOURS. 
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(b) LOWER WALL. 

FIGURE 46. - WALL PRESSURE DISTRIBUTIONS, CASE 2. 
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